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Abstract 

We present the perturbative solution of the multicomponent Boltzmann kinetic equation based 
on the set of observables including the hydrodynamic velocity and temperature for each component. 
The solution is obtained by modifying the formal density scaling scheme by Enskog, such that the 
density of each component is scaled independently. As a result we obtain the species momentum 
and energy balance equations with the source terms describing the transfer of corresponding quan- 
tities between different components. In the zero order approximation those are the Euler equations 
with the momentum and heat diffusion included in the form of the classical Maxwell-Stefan dif- 
fusion terms. The first order approximation results in equations of a Navier-Stokes type with the 
partial viscosity and heat conductivity including only the correlations of the particles of the same 
component. The first order corrections to the Maxwell-Stefan terms as well as the contributions bi- 
linear in gradients and differences of observables are calculated. The first order momentum source 
term is shown to include thermal diffusion. The nondiagonal (in component indexes) components 
of viscosity and heat conductivity appear as second order contributions. 

PACS numbers: 
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I. INTRODUCTION 



The theory of transport phenomena in gases based on a kinetic theory dates back to the 
works of Chapman and Enskog and is generally accepted as mature. It provides the 

perturbative solution of Boltzmann kinetic equation for the single- as well as multicompo- 
nent gases, resulting apart from other things in the microscopic basis for the Navier-Stokes 
equations. In spite of that the applicability of Chapman- Enskog solution to some particular 
problems is limited and a number of works has been published on modified and alternative 



derivation of transport equations from the microscopic kinetics Dynamics of 

multicomponent flows has also attracted considerable interest of the statistical mechanical 
community in the past. The works on the subject include the derivations of species mo- 



mentum and heat balance equations [8 
modification of diffusion equations JlOl 



transport theories can be found in 



, l9| , as well as the works devoted to discussion and 
Q, [13 1 . A recent overview of multicomponent 

ll 



The common recipe in modeling n-component transport problems is to solve the Navier- 
Stokes equation for the mixture, combined with {n — 1) Maxwell- Stefan equations, all with 
appropriate boundary conditions. Application of this principle to binary counterdiffusion 
in capillaries [l^ leads however to large discrepancies with experimental results (see for 
example 3])- On the other hand, instead of the set of observables including the number 
density of each component (nj), the hydrodynamic velocity (vq) and the temperature (T) of 
the mixture one may wish to use the set including the velocities and temperatures for every 



component 
and Stefan 



t^jn, Tj). This choice of observables was made already in the works of Maxwell 
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181], but this has become obscured by the success of the Chapman-Enskog 



theory. The species momentum and heat balance equations resulting from such an approach 
can be also considered as a basis for multifiuid models used in plasma theory 1^ . It should 
also be noted that that the tools to solve multifiuid transport equations numerically are well 
established 

Multifiuid models with application to the disperate-mass binary mixtures were studied 
before in the framework of kinetic theory. The two-fluid theory was developed by Goldman 
and Sirovich 2l| for the special case of Maxwell molecules, and later attempts were made to 



extend it to arbitrary interaction potentials 



22 1 . However, quite elegant but rather ad-hoc 



derivation of de la Mora and Fernandez-Feria 22] is restricted to binary mixtures, and, in 
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spite of the similarity of the resulting transport equations with the present work, the use of 
"unmatched Maxwellians" that are not a solution of any equation creates a lot of confusion 
throughout the treatment, especially when dealing with momentum and energy transfer be- 
;ween the components. The problem was recently approached also by Kerkhof and Geboers 



. Relating each component distribution function to its own hydrodynamic velocity, they 
obtained the momentum balance in the form similar to the Maxwell- Stefan equations. How- 
ever, in their derivation inaccuracies are present 23|] and thus serious reconsideration of the 
result is required. We like to stress that the present derivation is not related to ref. j?], al- 
though Kerkhof has also contributed to this work. We should also mention the little-known 



paper of Struminskii j2J| that we came across at the latest stages of development. To derive 
the equations for multicomponent gas transport Struminskii employed the same basic idea 
of modified scaling that is used in the present work. However our treatment is very different 
from given there, and we claim our results to be more general and accurate compared to 
those of Struminskii. More detailed comparison with this and other papers is given in the 
text where appropriate. 

In the present paper we aim to obtain the perturbative solution of the Chapman-Enskog 
type for the Boltzmann equation based on extended set of observables (nj, ViQ,Ti), and to 
establish equations of motion for these variables, relating the transport coefficients to the 
microscopic properties of the mixture components. The paper is organized as follows: In 
section II we suggest a modification to the standard density scaling scheme of solving the 
Boltzmann equation. Equations for the distribution functions as well as the momentum 
and energy balance equations are given there for the arbitrary order in the formal scaling 
parameter. In section HI we derive explicit expressions in the 0-th order in formal scaling 
parameter, while section IV contains the results of the higher orders of perturbation theory. 
Discussion of the results is given in section V. 

II. MODIFIED DENSITY SCALING 

We consider the multicomponent Boltzmann equation 

2)/. = J]J{M}, (1) 
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which describes the evolution T)fi = + Vr + Fj ■ 'V^Jfi of the component i distribution 
function /j due to colhsions with particles of the same and other components. The pairwise 
collisional operator is given hj J{fifj} = f dvjdk[j \vi-Vj\aij{\vi-Vj\,kij\k-j){f'J'j-fifj), 
where kij is the unit vector in the direction of relative velocity, and primed quantities denote 
those after collision. Here we only consider the case of spherically symmetric interactions, 
such that the crossection aij{\vi — Vj\, kij\k^j) to change the direction of relative velocity 
from kij to k'-j depends only on k'-^ ■ kij. We look for the solution of equation ([1]) based 
on the set of observables jSi = {rii, Vio, Ti) including the hydrodynamic velocities (■yjo) and 
temperatures (Tj) for each of the components. 

Following the standard procedure of Enskog we introduce the formally small scaling 
parameter e, put equal to 1 at the end of derivation. The solution is then given in terms of 
a series in this parameter /j = Xln^"//"^- -'■^ classical treatment this parameter is used 
to scale the left-hand part of equation ([T]) which is equivalent to scaling the density of the 
mixture. Such a scaling is very natural if one is going to use the set of observables including 
the hydrodynamic velocity and the temperature of the mixture as a whole (n^, Vq, T). Here, 
instead, we suggest a bit different scaling of equation ([1]) 

e[D/.-X;W.}] = W.}- (2) 

This corresponds to scaling only the chosen i-ih component density, and allows the zero 
order solutions to be centered around different velocities, which is convenient for deriving 
the species momentum and energy balance equations. It should be noted that the similar 



reordering was proposed yet by Lorentz 25| for the heavy species in disparate mass binary 



271, 



28| . However they 



mixture, and in modified form used later by many researchers [26, 
use such a reordering of the terms to collect all the relevant contributions in equations of 
the same order, thus requiring different scaling patterns for different components, and as 
a consequence, the different transport equations. Here instead, we use the same form of 
scaling for every component, employing it as a formal tool, which allows to construct the 
contributions from observable derivatives up to desired order. Namely, in the present paper 
we limit ourselves to the thermal conductivity, thermal diffusion and viscous contributions. 
Using the same form of scaling independently of component masses allows to treat the 
mixtures with arbitrary number of components, but on the other hand it leads to slower 
convergence of the series in rrii / Hij for some transport coefficients in the case of very light 



particles. 

Expanding unknown functions fi, operators D and J{fifj} as well as the time derivatives 
of observables df3i/dt in series in scaling parameter e, and assuming after Enskog that the 
time enters only implicitly through f3 and the spatial gradients of (3 



/,(r, V,, t) = /,(r, vf, (3, V,/3, . . .) = e"/, 



(n) 



d_ 
di 



A(r,t) = $.(r;/3,V,/3,---) = E^"^ 



(3) 



we arrive at the following equations for the contribution of the order n + 1 in a scaling 
parameter 

j^i k=i (4) 

^{/r/r } = 0, 

where the following notation is introduced 

n 

(5) 

^^"HM} = E'^{/^^/r"^}- 

m=0 

The right-hand side of equations (jlj) contains only the functions which are known from 
the lower order calculations. The left-hand side, containing the functions to be calculated, 
includes only coUisional operators for the species of the same type. 

The solubility condition for equations (jlj) , demanding that the right part of equation must 
be orthogonal to all solutions of the homogeneous equation, read as (the scaling parameter 
e is already put equal to 1 here) 

E / d'vii^.iv, -Vr + F- V,)/(") + E E / = 0. (6) 

n n j^i 

As long as the solution of the homogeneous equation (jlj) consists of collisional invariants 
ipi = {l,miVi, ^rriicf), with Ci = Vi — Vio the velocity related to coordinate set moving 
with the average i-th component velocity, the solubility condition provides us with the 



dt 
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corresponding conservation laws. After some manipulation this can be rewritten as 
1 dpi 

— -TT = -V ■ '"io. 
pi at 



n n n,j 



with partial pressure and the heat flow defined as 

(n) rrii 



(8) 



Here it should be noted that jniicf can be considered as a collisional invariant only in a 
limited sense as it conserves only under collisions of the particles of the same type. This leads 
to the sum of the sources I d^^i^^^^ J^'^^ {fifj} in the heat balance equation of the set ([7]) 
being not equal to zero. However there is no contradiction here with the energy conservation 
as the temperatures for each component are defined in its own reference frame, and thus the 
energy balance would have also included the terms of the pjf|)/2 type. Additionally, these 
source terms can only be calculated exactly in some special cases, and in what follows we 
expand them assuming the velocity Vio — Vjo and temperature Tj — Tj differences are small 
and keep only the first nonvanishing term of expansion. Thus, resulting equations should be 
understood as a linear in above mentioned differences approximation, while the sum of the 
source terms in (JTj) is obviously a contribution of the next order. Equations (171) can be seen 
as the single- component balance equations with the momentum and energy transport due to 
interaction between the components included as a source terms. This approach is equivalent 
to considering the single component system in a field of hydrodynamic forces exerted by the 
flow through the other components. We like to note however that quantities Pi and Qi can 

(n) 

depend on the properties of the other components through the functions f- . 

Similarly to original Enskog approach we derive here the transport equations in a general 
form with all contributions determined from equations (5,8). Calculation of these contribu- 
tions is then straightforward (although elaborate), with all assumptions and approximations 
being very clear. This should be contrasted to the previous derivations, which often require 
cumbersome argumentation to justify some steps in there. 
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III. ZERO-ORDER APPROXIMATION 



The solution of the equation (jl]) in the order approximation and resulting in the proper 
observables values is the Maxwell distribution related to the frame moving with the average 
velocity of the i-th component ■Ujo 

= n,(7./27r)3/2exp {-j,{v, - v,o)V2} , (9) 

where is the number density of the i-th component, 7^ = rrii/Ti with rrii being the mass 
of the component's molecules and Tj its kinetic temperature (we prefer using the energetic 
temperature scale, thus saving on ks). The resulting partial pressure and the heat flow are 



(10) 



where / is the unit tensor. 

The collision integral J^^^{fifj} can be expanded in series in velocity and temperature 
difference between the components. Then, assuming these differences are small, compared 
to the thermal velocity and the temperature of the component correspondingly, we keep only 
the first nonvanishing terms of expansion. The contributions of the zero order in velocity 
and temperature differences are identically zero, while the most general form of the first 
order contribution would be 

J^'^im = -fPQ^'Kcdiv^o - v,o) ■ c. - /f g(^)(Q)(T, - T,). (11) 

Explicit form of the scalar functions Q^'"\ci) and Q^'^^Ci) is rather complicated, and for the 
time being we keep them undefined. In the end it turns out that we only need to know 
the even moments of these functions with respect to q and not the functions themselves. 
Calculation of these moments, using equation ( ITTi) gives 

f ,3^. 2n.{0)^(T). ^ _ _^r(2n + 2) /^^X"! ^ 7. + ^Ij'"^ 



n n+1 (1;™.) 



J ' Jl ^ \i) *22"+in+l V 2 ^ 

\ / \ / m — 1 



r(n — m + 2) ' 

(12) 



where the hierarchy of frequencies is introduced 



„ ™02m+3 /,, \ 2m-l / \ m-l ^oo 



with ^ = 7jj(yf/2, fiij is the reduced mass of the components i and j, 'jij is similarly introduced 



reduced quantity, and Qfj{g) = j dQaij{g,6){l — cos' 6*) is the ij transport cross section of 
the l-th order. Q'l^\g) entering the zero order collision integral is sometimes referred to as 
a diffusion or momentum transfer cross section. When all components of the mixture have 
the same kinetic temperature these frequencies can be easily related to familiar Chapman- 
Cowling fi-integrals: u^j'"^^ \Ti=Tj = (rrij /mj)"^( tijj / m,- ) ^2^''"^'* . The finite series for the 



moments of collision integrals were calculated by Kolodner 



29| for the special case of Maxwell 



molecules, and were used by Goldman and Sirovich j2l| to develop the two-fluid theory for 
the mixtures of those molecules. Here we use the similar formalism to derive the multifiuid 
theory for the species with arbitrary interaction potential. This should be contrasted to 



the work of Struminskii 2J], where similar quantities were also calculated for the case of 
Maxwell molecules, and implicit assumption made that the same expression can be applied 
to the species with arbitrary interaction potential. 

Extracting the required moments of J^^^{fifj} from f fT2l) we can write down the transport 
equations in the following form 
1 dp^ 

— -1- = - V • Vio, 
Pi at 

dvio (i-i) 
Pi^jf = PiFi - ^Pi - ni fiijUlj '{Vio - Vjo), ^^4^ 

j 



The I'i^'^^ frequency here coincides with the effective coUisional frequency of the electron 



commonly used in the plasma theory 19|. Thus in the zero order approximation we obtain 
the Euler equations for the ideal liquid including the Maxwell-Stefan terms for the heat 
and momentum diffusion due to interaction between different components. The zero order 



results are in total agreement with Struminskii 24| and binary Mora 22(] equations, allowing 
to recover the relationship between diffusion rates, obtained by the Enskog method in the 
first approximation. 

IV. FIRST ORDER APPROXIMATION 

The first order contribution to the distribution function is conveniently put in the form 
/j*-^-* = ff^4>t\ introducing the new unknown function ^-^^ Using the results of the previous 
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section to calculate and substituting expression for J^^^{fifj} we can rewrite equation 

(jlj) for the unknown function (pi as follows 



(0) 



1 



It 2 



-cf n : V^;,o +\^ci--\ c, ■ VlnT, 



3 



mi + mj ^ 



where I{(f)} = J d'^cidQ\c — ci\f^^''fi^\(j) + (f)i — (j)' — (j)[) is the linearized collision operator 
[3]. This expression differs from the similar single component equation only by the 'friction' 
terms proportional to {ViQ — Vjq) and (Tj — Tj), therefore the corresponding contributions 
should be included into the function (h^^\ The structure of this expression is similar to that 

n □ 

used by Goldman and Sirovich [211], however in the later work of de la Mora [22| the Q-terms 
arising from collision integrals were omitted. Taking into account rotational invariance of 
collision operator we can write the function (pi in the form 

0i = -—A{ci)Ci ■ VlnTi - —B{ci) ( CiCi - ^c^J ) : Vvio 

Hi rii \ 3 J 

- A,(q)c, ■ {v,o - v,o) -- Y,E,,{cim-T,). (16) 



rii '-r. rii ... 



The form of the function is similar to that used in the work of de la Mora but differs 



from the one in the paper of Struminskii 



3, 



2J] , where the unknown functions were related not 



to the driving forces, or independent parameters of solution, but rather chosen in accordance 
with tensorial order of Cj weights. In our opinion this can not be done if one is looking for 
solution based on extended set of observables. 

In order for (p^p to satisfy equation f[T^ the following identities should hold 

nd{A{ci)c,} = f^- 



.2 5 
2 2 



(17) 



niI{B{ci){c,c, - ^cU)} = l^f'\ciC, - Ic^/), 

.wcO}^/^°)(q^-(.) + |0-^^ 

Since the proper values of all observables are obtained from the 0-order distribution func- 
tion, several conditions should be apphed to the function 0j, such that f^^^ provides no 



contribution to observables. In terms of functions A{ci),B{ci) and E{ci) these conditions 
read as 



Jd'cc'f^'^D,,ic) = 0, (18) 

which determines the function 0j uniquely. 

To determine the functions A{c), B{c), D{c), E{c) we use the variational approach [s], Q, 
i.e. we look for the functions that maximize the corresponding bracket integrals and satisfy 
(I17fl8l) . Here we should mention that unlike in the classical treatment it is not immidiately 
obvious that the entropy production would also be maximized by this procedure, and the 
links between the current development and nonequilibrium thermodynamics require a sep- 
arate study. Expanding unknown functions into series we can obtain the simple equations 
to determine the expansion coefficients. The functions A{c) and B{c) here are exactly the 
same as in classical single-component approach, thus resulting in the 'partial' viscosity and 
thermal conductivity, equal to those the i component would have by its own at the same 
number density rii and temperature Tj it has in a mixture. Expressions for these functions 
are well known, but we provide them here again for completeness. Similarly to A{c) and 
B{c) the functions D{c) and E{c) are most easily determined if we expand them in series of 



associated Laguerre polynomials [3|, l3ll] 



(19) 



A,(c.) = -(7./2)V2 5^a;;)43/2)^^^^2/2), 

n=l 

oo 

S.(c.) = (7./2)5:6y4^/^n7AV2), 

n=0 

oo 

A,(q) = -(7./2)5^<)L(f/^)(7.c?/2), 

n=l 

oo 
n=2 

Some terms in these series are omitted here to satisfy conditions (|T8l) . 

Using the moments of the the functions Q^'"\ci) and Q^^\ci) ( fT2l) we can write down the 
simple matrix equations to determine the expansion coefficients e[^^ and d^™^ (the equations 
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4. 



for Oj and bi coefficients are again provided here for completeness) 

oo 19 

J2 r^[(7./2)(c.c, - -c?/)L(^/^n7.c?/2), (7./2)(c.c, - -cU)L^:/'\7^cl /2)] = -5„o 

m=0 * 

oo 



m=l 



oo -. 

(20) 

where the bracket integral is defined as [-F, G\ = Tr J d?Ci G ® I{F} and the shorthand 
notation K^f = ELo (l - ^)"" (^)' - introduced. 

Provided equations (|20l) are solved and thus the functions A{c), B{c), D{c) and E{c) are 
known, we can calculate the first order contributions to equations (171). For the pressure and 
the heat flow those read as 



5 (1)^^ 5T,^ ,(1), , (21) 

— ( 
4 



Here, the partial viscosity and heat conductivity are given by the ffist expansion coefficients 
bf^Ti/2 and ^af'^ /A correspondingly. There is also a term describing the heat flow due to 
interaction of the flows moving with different velocities, however temperature differences do 
not contribute to the pressure. Instead there appears a term in the ffist order momentum 
source integral having the same effect. The partial viscosity and heat conductivity here are 
determined by only the ffist nonvanishing coefficients from the expansion (fT9|) . On contrary 
the contributions from the friction terms, are given in terms of series including the whole 
list of expansion coefficients. Up to the terms bilinear in gradients of observables and their 
differences the source terms in the ffist order calculate 
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and 

-Y.^>^i',liT^-T,) + \f,{T,-T,)). (23) 

k 

The transport coefficients /i*^) and A^^ include the whole expansions {a^")}, {ft*-"-*}, {(i*^")} and 
{e(")}, and are given in appendix A. The bilinear terms here are not really important and 
most probably can be omitted for most of the practical applications. The most interesting 
contribution is probably the one describing the thermal diffusion which explicitly appears 
as a result of interactions between different components in the current treatment. If the 
components share the same value of mj and Tj (that is for example if the single-component 
flow is virtually divided in two flows of similar species), this contribution disappears in total 
agreement with the classical solution. Summing the equations to obtain the momentum 
balance for the mixture as a whole also eliminates this term. 

The other interesting contributions are those with transport coefficients with 3 indexes. 
Those not only include the ffist order corrections to the Maxwell- Stefan diffusion terms, 
but also describe indirect interactions of two components through the impact on a third 
one. Appearance of such contributions derived from the Boltzmann equation, allowing only 
the pairwise interactions is fascinating but not totally unexpected. Use of the reduced 
description with some degrees of freedom integrated out often leads to many-body effective 
interactions. The most famous example of that is probably the interaction of the ions of the 
same sort in electrolytes. 

Thus, in the ffist order approximation we obtain the Navier-Stokes equations which in- 
clude the corrected Maxwell-Stefan terms, thermal diffusion and also several new terms bi- 
linear in gradients and differences of observables. The partial viscosity and heat conductivity 
depend only on the properties of the component thus totally neglecting ij correlations. 
Also the time evolution of Vio depends on the strain tensor of the component i and does 
not include the strain of the other components, which is not physical. This situation can be 
understood if we notice that the formal density scaling ([21) we are using, leads to mixing the 
terms of different orders of the scaling parameter e as compared to the standard treatment, 
and tends to give more weight to ii interactions. Thus, although, the second order distribu- 
tion function does not contribute to viscous forces in the classical approach, here we need 
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to proceed to higher orders to include the ij viscous interactions. 

As many contributions from the first order treatment can aheady be considered excessive, 
we are not going to do the complete analysis of the second order solution. Instead we 
are only interested in correcting the viscosity and heat conductivity including the terms 

(2) 

resulting from ij correlations. As previously, introducing the auxiliary function 0^ such 
that /j = f- we rewrite equation (jlj) in the form 



dt dt 



Here Ylj^ti •^^^H/j/i) is the only term containing the required contributions, the most general 
form of which is given by 



(0) 



(25) 



+Q^;^^(q)(qc, - icf/) : S, + Q^J'^'ic^c, . VT, 

We use the sign oc to indicate that all irrelevant contributions are omitted here. The func- 
tions Q^jj{ci) can be left undefined for the time being, as for the later analysis we will only 
need their moments with respect to even powers of q which are easier to calculate then the 
functions themselves. 

Similarly to the previous section, the viscous contribution to the second order pressure 
tensor is found 

i^^^^ = -^EC^-^E^r^' (26) 

and nondiagonal part of the heat conductivity reads as 

<f>--lE4''^^-lj:^''^T, (27) 

The expansion coefficients again can be found from the simple matrix equations 



oo 

Y: ag-)[(7./2)V2c.L(;^/^)(7.c,V2), (7./2)V2QL(f/^)(7.c?/2)] = K 



(28) 



Expressions for the coefficients K^j'""^ are given in Appendix B. When the mixture is treated 
as a whole and thus all components have approximately the same rate of shear and temper- 
ature gradient, the viscosity and heat conductivity for the mixture can be found simply as 
a sum of all the components of corresponding matrices. 
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V. DISCUSSION 



We investigate the perturbative solution to the Boltzmann multicomponent kinetic equa- 
tion based on the set of observables including each component velocity and the temperature. 
The corresponding set of the species balance equations is derived. The new solution is ob- 
tained through modification of the formal density scaling scheme of Enskog, such that the 
density of each component is scaled independently. This approach allows for succesive devel- 
opment, with clear assumptions and approximations, which should be contrasted to rather 
ad-hoc previous attacks on the problem. Omittance or addition of the terms in some inter- 
mediate equations as compared to the similar works, as discussed in the text, follows directly 
from the scaling equation 1^. 

We show that the distribution functions related to each component velocity appear nat- 
urally using this approach, and the solubility conditions immediately provide us with the 
balance equations for species. The zero order equations are the Euler equations for each 
component including the heat and momentum transport between components through the 
familiar Maxwell- Stefan diffusion term. As a next approximation we obtain the Navier- 
Stokes equations with partial viscosities equal to those the component would have in ab- 
sence of the other. Apart from that there are viscous contributions from the shear of other 
components. The interaction between components leads to appearance of the heat diffusion 
and also some corrections to the Maxwell- Stefan diffusion term. These corrections include 
apart from direct interactions between the two components indirect ones, through perturb- 
ing the other components flows. The heat balance equations posses the similar properties. 
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To summarize we provide here the final equations describing the multicomponent fiow 
1 dpi 

--TT = - V - ■Uio, 
Pi at 

j k 

= - P»V ■ Vio + VA.VT, + 5^(VAi,- VT, + VA,-,VT,) 
+ 2?7i5i : VVio + 2 ^(r/jj-Sj + rijiSj) : Vvjo 

(29) 

where the transport coefficients are deifned in the text through the propeties of the com- 
poenents. 

The convergence of the scheme is not questioned in the paper, however one may argue that 
apart from using an extended set of observables the whole procedure can be seen as mixing 
the terms of different orders (compared to the classical treatment), and thus convergence of 
both solutions should be about the same, suggesting the range of Knudsen numbers where it 
is applicable. Another issue is the convergence of the transport coefficients which depend on 
the whole first order distribution function, and not just include the first expansion coefficient. 
As we see, convergence of these series depend strongly on the value of < ■Jij/'jj < 1. In 
cases when the temperatures of the different components are about the same, and the impact 
of the temperature gradient is more significant, which is true for many chemical engineering 
applications, the above ratio can be replaced by pij/nij, thus removing the complicated 
temperature dependency from the transport coefficients. In the case when the gradients are 
small, but the temperatures themselves are very different, as it may happen for example 
in plasma, the Jij/'jj ratio should be kept, but it varies slowly again. The only potential 
problems may appear when both gradients and temperature differences are big enough, 
leading to complicated temperature dependency of the transport coefficients, however this 
case is clearly out of scope of the present theory. The present study can also be seen as 
a basis for the multiffuid model used to describe kinetic processes in plasma. The other 
potential applications of the method presented here include reactive ffows and the boundary 
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layer problem for mixtures o 



literature (see for example 



' rarified gases, that still generates a lot of discussion in the 



32 



33|). Some authors also tend to employ the multicomponent 



transport equations to describe the shock- waves or even turbulence phenomena 34| |. 

Apart from modified scaling, the presented solution follows closely the original Enskog 
scheme, with which it has more in common then any other work on multicomponent trans- 
port cited here. However some features appear to be very different. The question of entropy 
production which is not obviously maximized by the present solution requires separate study 
to establish the links with the nonequilibrium thermodynamics. The fact that there are no 
ij collisions in the left part of the equation (jl]) eliminates the need of ij bracket integrals. In 
some sense they are replaced by the set of collisional frequencies z/^™"'"^ , closely related to the 
Chapman-Enskog ^-integrals. The J^lj'^^ is recognized as an effective collisional frequency of 
electrons from the plasma theory Unlike in the classical solution the first order solution 
does not provide us with all the relevant contributions, and the contributions of the same 
nature, such as viscous, are divided between the first and the second order contributions. 
This is the consequence of the suggested scaling, which gives more weight to collisions of 
the species of the same type. On the other hand the ii and ij contributions to the viscosity 
and heat conductivity are explicitly separated, which makes it convenient to compare the 
present results with the Green-Kubo theory for the transport coefficients. It should also 
be noted that for the special case of binary mixtures the form of the transport equations 
derived here coincide with those obtained by de la Mora 22] . The structure of the equations 
for the transport coefficients however is somewhat different, although we do not expect large 
discrepancies in the numerical values. The detailed numerical study revealing importance of 
different contributions, and containing comparison with counterflow and other experiments 
(simulations) is the topic of the future research. 
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APPENDIX A 



Below we provide the expressions for the transport coefficients entering the first order 
equations of motion. Lower signs correspond to switching i and j indexes, and F{a, b; c; x) 
is the hypergeometric function. 
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APPENDIX B 



Here we provide the expressions for coefficients appearing in the right hand side of equa- 
tions (I28p which determine the nondiagonal (in component indexes) viscosity and heat con- 
ductivity. Here again the ij coefficients are obtained by simultaneous replacing all i by j 
and contrary, no sign changes required here. These expressions are rather awkward, and 
in principle several sums should be possible to calculate analytically, however we failed to 
further simplify these expressions, which anyway would not make significant difference for 
the numerical calculations of the coefficients. 
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Here we have introduced an auxiliary function Wf"^ „ 



also implied that only contributions with integer upper indexes of u^j enter the summation. 
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